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Abstract. We present CO 1— >0 and 2— »1 observations of 
the dwarf starburst galaxy NGC 1569 with the IRAM in- 
terferometer on Plateau de Bure. We find the CO emission 
is not spatially associated with the two super star clusters 
in the galaxy, but rather is found in the vicinity of an HII 
region. With the resolution of our data, we can resolve the 
CO emission into five distinct giant molecular clouds, four 
are detected at both transitions. In the 1—^0 transition the 
sizes and linewidths are similar to those of GMCs in the 
Milky Way Galaxy and other nearby systems, with diam- 
eters ranging from ~ 40 to 50 pc and linewidths from 4 
to 9 kms" 1 . The (2-l)/(l-0) line ratios range from 0.64 ± 
0.30 to 1.31 ± 0.60 in the different clouds. The lower line 
ratios are similar to those seen in typical Galactic GMCs, 
while values higher than unity are often seen in interacting 
or starburst galaxies. We use the virial theorem to derive 
the CO-H 2 conversion factor for three of the clouds, and 
we adopt an average value of 6.6 ±1.5 times the Galactic 
conversion factor for NGC 1569 in general. We discuss the 
role of the molecular gas in NGC 1569, and its relationship 
to the hot component of the ISM. Finally, we compare our 
observations with blue compact dwarf galaxies which have 
been mapped in CO. 

Key words: ISM: molecules - Galaxies: individual: 
NGC1569 - Galaxies: ISM - Galaxies: starburst - Radio 
lines: galaxies 



1. Introduction 

NGC1569 (Arp 210, VII Zw 16, UGC3056) is a nearby 
dwarf galaxy hosting several interesting phenomena re- 
lated to its starburst, and with an observational history 



going back to 1789 (see Israel 1988 for a history of the 
early observations). As is common among dwarf galaxies, 
it has a low metallicity (12 + log(0/H) = 8.19 ± 0.02; 
Kobulnicky & Skillman 1997). At a distance of only 2.2 



± 0.6 Mpc (Israel 1988), it is the closest known example 



essential in interpreting observations of similar objects at 
greater distances (e.g. the small blue galaxies found in the 
Hubble Deep Fields). 

What has perhaps drawn the greatest attention to 
NGC1569 is the presence of two super star clusters (SSCs), 
labeled A and B (Abies [l97T| , Arp & Sandage |l985| ) . These 
clusters have been th e subj ect of recent HST ( O'Co nnell, 
Gallagher & Hunter |l994[ De Marchi et al. |1997|) and 
groun d based studies (Prada et al. 1994 ; Ho & Filippenko 
1996 ; Gonzalez Delgado et al. 1997 ), and are believed to be 
similar to young globular clusters. Age estimates for the 
SSCs range from 3 to 10 Myr depending upon assumptions 
about the star formation history. 

Greggio et al. ( |1998| ) have determined from HST 
WFPC2 images that NGC 1569 has experienced a global 
burst of star formation 100 Myr in duration that has ended 
as recently as 5 Myr ago. Vallenari & Bomans (1996) 



found evidence for a large burst of star formation roughly 
1 x 10 8 yr ago in WFPC images, as well as several much 
older episodes. It is clear that star formation has had a 
dramatic effect upon the ISM in NGC 1569 through stel- 



lar winds and supernovae. Israel & van Driel (1990) have 



found a hole in the HI distribution centered on SSC A, 
possibly blown out by the stars in the cluster. Ha emission 
extends from the di sk out t o the halo in filamentary struc- 
tures ( Hodg e |1974| ; Waller |l99l[ H unter Hawley & Gal- 
lagher |1993|; Tomita, Ohta & Saito [1994 Devost, Roy & 
Drisscn 1997 ). The dynamical age of the extended, diffuse 
Ha emission is consistent with some age determinations of 
the S SCs (H eckman et al. 1995). X -ray studies (Heckman 
et al. 1995 ; Delia Ceca et al. 1996 ) have found extended 
emission spatially associated with the Ha; together these 
paint a picture of a hot gaseous phase blowing out of the 
galaxy, powered by the SSCs. 

The cool phase of the ISM in NGC 1569 has been ob- 



of a dwarf starburst galaxy, and so observations of it are 
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served by Israel & van Driel ( |1990| ), Stil & Israel fll998| ) 
and Wilcots et al. (in preparation) in HI with interferom- 
eters. Hunter et al. ( |1989| ) have combined FIR data from 
the Kuiper Airborne Observatory and IRAS to study the 
dust. They find an unusually high dust temperature, for 
a dwarf irregular galaxy, 34 K, and attribute this to the 
influence of the recent star formation burst. CO emission 
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has been observed in NGC1569 by Young, Gallagher & 
Hunter Jl984| ), Greyeet al. fll996| ) and Taylor, Kobul- 
nicky & Skillman ( |1998[ ). Mapping with the IRAM 30-m 
telescope, Greve et al. found CO emission near the SSCs, 
but not directly at their locations. Instead the CO is spa- 
tially associated with a prominent HII region. The spatial 
resolution of the 30-m telescope is 22" at 115 GHz (equiv- 
alent to 235 pc) and 11" at 230 GHz (equivalent to 120 
pc), so it was not possible to resolve individual molecu- 
lar clouds. Because of this, an extremely high value of the 
CO-to-H2 conversion rate, ~ 20 times the Galactic value 
was derived. Such a high value is unlikely, and Greve et al. 
attributed it to non-standard conditions in the molecular 
gas due to the proximity of the super star clusters. 

NGC 1569 presents an excellent opportunity to study 
the effects of low metallicity and a starburst upon the 
molecular ISM. The only similar environment which can 
be studied in such detail is the molecular gas complex 
found near 30 Doradus in the LMC. We have imaged the 
CO emission using the IRAM interferometer on Plateau de 
Bure, which provides spatial resolution sufficient to distin- 
guish individual giant molecular clouds (GMCs). We ex- 
amine the physical properties of the GMCs, derive the CO- 
to-H2 conversion rate, and discuss the role of the molec- 
ular ISM in NGC 1569 with respect to its star formation 
history. 



64 51 15 

00 

50 45 



NGC1569:CO(1-0) 

-61,34<iM/S 



64 51 15 



> o_ 
oo o ° 



50 45 






—64 51 15 

z 
o 

< 00 



o 

g 50 45 



64 51 15 

00 

50 45 



-:0' : 



I- 



o 






<?<©' 



-# 



■ti.Vkm/s 



0-< 



o 



<g0.8 KM/S 

6° 










-90i6 KM/S 



O 



-64.6K1WS 



.0^ 






«7*3 KM/S 

^ " o Q 



,®>„a 



JrtKM/s 



h — "y " ., e '■: t 



V93M KM/S 

f ° '»V 



3 






t~ 






-&7.8Kjg/S 



; m "' s 



~4- 



.-V77.6KB/S 

"o •■■■'; 



v 



187.9|E)*S 



.•i "s? ) °c y y^ 






c-p '.; n _ -1 



* r o t97:1 K^/S 



^?"! 



a 



,e>- 



04 30 50 48 46 44 

RIGHT ASCENSION (J2000) 



04 30 50 48 46 44 



2. Observations and Data Reduction 



The observations were carried out with the interferometer 
on Plateau de Bure in several sessions between 21 March 
and 14 April 1998. The 12 CO 1-vO and 2->l transitions 
were observed simultaneously in a mosaic consisting of 
four positions separated by 11", centered upon the peak 
in the CO emission detected by Greve et al. ( 1996J ). The 
D and C2 configurations were used, yielding spatial reso- 
lutions of 4"5 x 3"9 at 115 GHz and 2"3 x 2"0 at 230 
GHz. At the adopted distance of 2.2 Mpc, these spatial 
resolutions correspond to 45 pc x 40 pc and 25 pc x 20 
pc, respectively The velocity resolutions are 1.6 kms -1 at 
115 GHz and 3.3 kms" 1 at 230 GHz. 

The data were reduced at IRAM with the CLIC and 
MAPPING packages of GILD AS, using the standard pro- 
cedures. During the observations the atmospheric water 
vapor content varied, and was often high enough to influ- 
ence the 230 GHz observations. The phase calibration was 
accomplished accurately with the help of the stable phases 
at 115 GHz, but the amplitude calibration at 230 GHz is 
only accurate to 40%. The accuracy of the 115 GHz am- 
plitude calibration is 20%. After mapping was completed 
the data cubes were exported as FITS files to the AIPS 
package for further analysis. 

The 115 GHz data have an rms noise of 7.5 mJy/beam 
in a single channel map, while for the 230 GHz data the 
rms noise is 6.6 mJy/beam. The channel maps are pre- 
sented in Figures 1 and 2. 



Fig. 1. Channel maps showing 12 CO 1—^0 emission in 
NGC 1569 near SSC A and B. The contours are in units 
of -2, 2, 4, and 8<r, where a is the rms noise in a channel 
map, equal to 7.5 mJy beam -1 . Every second channel is 
shown. Velocities are vlsr- 



Each data cube was blanked at the 3a level and the re- 
sulting blanked data cubes were searched for CO emission. 
To distinguish genuine emission from noise spikes, addi- 
tional blanking was done in which only emission present in 
at least three consecutive channels was retained. The inte- 
grated CO emission over the entire map is 3.58 Jy kms -1 
in the 1— >0 line, and 5.94 Jy kms" 1 in the 2— >1 line. In 
comparison, Greve et al. ( 1996| ) integrated the emission 
over the inner 22" of their map, obtaining 12.5 Jy kms -1 
(1.98 K kms" 1 ) and 22.1 Jy kms" 1 (2.10 K kms" 1 ) for 
the 1— >0 line and 2— >1 line, respectively. Our field of view 
is larger than this, so if we restrict ourselves to the same 
inner 22 , we obtain 2.76 Jy kms and 5.79 Jy kms . 
At 115 GHz, we detect ~ 22% of the flux from Greve et 
al., and at 230 GHz we detect ~ 26%. This discrepancy 
is most likely due to the incomplete coverage of the uv 
plane at the shortest baselines. To check that insufficient 
sensitivity was not the reason, we assumed that the en- 
tire 22" region in our 1— >0 map, excluding only the areas 
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Fig. 2. Channel maps showing 12 CO 2— >1 emission in 
NGC 1569 near SSC A and B. The contours are in units 
of -2, 2, 4, and 8c, where er is the rms noise in a channel 
map, equal to 6.6 mJy beam -1 . Velocities are Vlsr- 



covered by the GMCs, was filled with emission 3 chan- 
nels wide at the 2.5<r level, i.e. just below our sensitivity 
defined above. If this were true, then the total emission 
should be 3.82 Jy kms -1 , and we would have recovered 
only 72% of this. But even this most extreme scenario still 
falls short of the 12.5 Jy kms -1 seen by Greve et al. Thus 
we conclude that the discrepancy is due to the missing 
short spacings. This lack of short spacings means that our 
observations do not detect diffuse gas which is distributed 
over large scale lengths, but only the dense gas which has 
accumulated into giant molecular clouds. Our results will 
therefore only be applicable to the structures we see, and 
not to the global distribution of CO in NGC1569. 

3. Results 

3.1. Cloud Diameters, Line Widths and Line Ratios 

The zeroth moment maps, showing the velocity-integrated 
emission, are given in Figures 3 and 4. Figure 5 shows the 
line profiles for the clouds in both CO 1— >0 and 2 — s-1 emis- 
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Fig. 3. Integrated CO intensity map of the 12 CO 1— >0 in 
NGC 1569. The contours represent 10, 20, 40 and 80% of 
the peak integrated flux intensity, equal to 1.6 Jy kms . 
The positions of the two SSCs are labeled A and B. The el- 
lipse indicates the size and shape of the synthesized beam. 



sion. Although the CO emission from some giant molec- 
ular clouds overlaps spatially, the velocity information 
shows that several distinct clouds may be distinguished, 
especially in the 2— >1 data, which have a higher spatial 
resolution. Four clouds are identified in each data cube, al- 
though clouds 1 and 2 from the 2— >1 are merged together 
by the low spatial resolution in the 1— >0 transition. Cloud 
5 from the 1— >0 data falls outside the area covered in the 
2— <fl map and thus is not seen. The observed properties 
of the GMCs are given in Table 1 , including the diameter 
measured at the contour encircling 90% of the flux (used 
for determining the virial masses) , the diameter measured 
at the half maximum contour (used for comparing with 
the size-linewidth relation), the FWHM velocity width, 
the central velocity, and the integrated intensity. 

One question to consider is how these GMCs in 
NGC1569 compare with those from the Milky Way 
Galaxy, and from other nearby galaxies. Both in the Milky 
Way and other galaxies, GMCs are observed to follow a 
size-linewidth relation of the form v oc D 13 where v is 
the linewidth, D the diameter, and /3 ~ 0.5 (e.g. Larson 



1981, Solomon et al. 1987, Wilson & Scoville 1990). Fig- 



ure 6 plots our clouds in the size-linewidth plane along 
with clouds observed in M31 (Vogel et al. 1987; Wilson 
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Table 1: Observed Properties of GMCs in NGC1569 



a(2000) 5(2000) 



D90 
(PC) 



■L' / w hm 
(PC) 



V fwhm 

(kins" 1 ; 



V center 

(kms" 1 ) 



S(CO) 
(Jy kms~ 



04 30 46.1 +64 50 55 110 x 57 70 x 44 

04 30 46.6 +64 50 58 72 x 57 61 x 35 

04 30 47.9 +64 51 08 75 x 45 55 x 33 

04 30 47.2 +64 51 14 81 x 56 55 x 37 



24.1 


-75.7 


2.12 


8.8 


-79.3 


0.86 


3.9 


-85.0 


0.19 


5.7 


-92.1 


0.41 


11.3 


-68.7 


1.98 


13.1 


-69.1 


2.28 


11.4 


-79.9 


0.83 


4.9 


-84.9 


0.23 



2- 

1 
2 
3 
4 



04 30 45.8 +64 50 54 41 x 34 31 x 24 

04 30 46.1 +64 50 55 47 x 38 37 x 28 

04 30 46.7 +64 50 59 49 x 26 36 x 21 

04 30 47.8 +64 51 08 28 x 21 24 x 16 



Note to Table 1: 

diameter ± 24 x 

CO 2^1 data 



The estimated errors on the various properties in the CO 1— >0 data are: 
21 pc, v fwhm ± 1.3 kms -1 , v center ± 1.3 kms" 1 , I(CO) ± 20%. For the 
the estimated errors are: diameter ± 13 x 11 pc, v/ m / lm ± 3.2 kms -1 , 

Vcenter ± 3.2 kfllS^, I(CO) ± 40%. 



NGC 1569: CO (2-1) 
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Fig. 4. Integrated CO intensity map of the 12 CO 2— >1 in 
NGC 1569. The contours represent 10, 20, 40 and 80% of 
the peak integrated flux intensity, equal to 1.8 Jy kms -1 . 
The positions of the two SSCs are labeled A and B. The el- 
lipse indicates the size and shape of the synthesized beam. 
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(Rubio et_aL |1992|), IC10 (Wilson |l995| ) and NGC 6822 



(Wilson 1994 ) . The li ne sho ws a fit to the points from M33 
by Wilson & Scoville ( |l990| ) of the form: v = 1.2D 05 . Our 
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Fig. 5. Line profiles for GMCs 1, 2, 3, 4 and 5. Clouds 1 
and 2 are blended together in the 12 CO 1— >0 data because 
of the lower spatial resolution. Cloud 5 is only shown in 
the 12 CO 1— >0 transition because it lies outside the area 
imaged at 230 GHz. 



clouds clearly fall within the range spanned by the clouds 
from other galaxies. 

The CO (2->l)/ (l->0 ) line ratio, r 2 i, is given in Ta- 
ble 2. Greve et al. ( 1996[ ) determined the line ratios for 
the inner 22" of their single dish map, obtaining 1.1 ± 
0.2. This value lies between the extremes that we have de- 
termined for the individual GMCs, although our error bars 
are large enough to include their value. If we integrate only 
the emission corresponding to the same region as Greve et 
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Fig. 6. The diameter-line width relationship for GMCs 
in galaxies of the Local Group. The clouds of different 
galaxies are distinguished by different symbols, and the 
line shows a best fit rel ations hip to the M33 GMCs derived 
by Wilson & Scoville fll990|). 



Table 2: CO (2-»l)/(l-»0) Line Ratios 



GMC 


2->l/l->0 


1+2 


1.31 ± 0.60 


3 


0.64 ± 0.30 


4 


0.79 ± 0.36 



al., we obtain a line ratio of 1.38, larger than their value, 
but consistent, given our errors. However, GMCs 1+2 con- 
tain most of the emission at both frequencies, so it is likely 
that the value of Greve et al. is dominated by the com- 
bined contribution of these two objects, and that the lower 
line ratios we obtain for GMCs 3 and 4 are correct. 

3.2. The CO-H2 Conversion Factor 

Under the assumption that the clouds we observe are in 
virial equilibrium, we may use the virial theorem to calcu- 
late the mass within the clouds. Since the mass of GMCs 
is dominated by molecular hydrogen, comparing the mass 
thus obtained from the observed CO flux density will give 
the conversion factor between CO and H 2 . For this calcu- 
lation to apply, the clouds must be resolved both spatially 
and in velocity, otherwise upper limits on the molecular 
mass of the GMCs are obtained, resulting in lower limits 
on the conversion factor. Figure 5 shows that the veloc- 
ity resolution is sufficient, although we will not be able to 



use the CO 1— >0 data for clouds 1 and 2 because they are 
blended together at the lower spatial resolution. At a dis- 
tance of 2.2 Mpc, the spatial resolution at 115 GHz of 4'.'5 
x 3'.'9 corresponds to 48.0 x 41.6 pc, while the resolution 
at 230 GHz of 2 '.'3 x 2 "0 is equal to 24.5 x 21.3 pc. 

An important question is whether or not the assump- 
tion of virialization is justified. From Figure 6 we conclude 
that the clouds are similar to GMCs known in the Galaxy 
and other nearby systems, so if those clouds are virialized, 
we may reasonably assume that ours are as well. That 
GMCs in the Milky Way are virialized has often been the 
subject of vigorous debate. One line of reasoning that ar- 
gues for virialized clouds has been the tight correlation 
between virial masses, M77, and CO luminosities, Lcoj 
where \jco is taken as an indicator of the mass of CO, and 



hence H2, present i n a gi ven cloud (Solomon et al. 1987 ). 
However, Maloney ( 1990 ) has argued that this Mvt~^co 
correlation is simply a result of the observed size-linewidth 
relation and would exist whether or not the GMCs were in 
virial equilibrium. In the end, the fact that conversion fac- 
tors derived with this method generally agree with those 
derived from independent methods, at least at the high 
mass end of the GMC distribution (> 10 5 Mq), suggests 
that the assumption of virialization is reasonable (Combes 
1991). Both sensitivity and spatial resolution limit our ob- 
servations to these largest, most massive GMCs, so we will 
use the viral masses to derive the H2-CO conversion fac- 
tor. 

The virial mass contained in a cloud is given by 



M VT = 190 



l! 



fwhm 

kms -1 



D/2 
pc 



AL 



(1) 



where v 2 fwhm is the velocity width and D the diameter of 
the cloud (MacLaren et al. 1988). For the diameter, we 



use the average of the major and minor axes, measured 
at the contour containing 90% of the flux (D90), and for 
the velocity width, we use the linewidths in the 1— >0 line. 
The factor of 190 is appropriate for a spherical distribution 
with density proportional to 1/r. The molecular mass in 
a given cloud of integrated CO 1— >0 flux density Sco [Jy 
kms _1 l 



is: 



M mol = 1.23 x 10 4 



-Sec 



Mpc Jykms" 



Mr. 



(2) 



wher e d is the distance to the cloud (Wilson & Scoville 
1990 ). This formula uses a Galactic conversion fac tor o f 
ct Ga i = 2.3 x 10 20 cm" 2 (K kms- 1 )- 1 (Strong et al. |l988|) 
and includes the helium correction. The CO-H2 conversion 
factor is then obtained from 



Al- 



oe = aGal 



VT 



M, 



raol 



(3) 



Table 3 lists the values for aj a ga i and Mvt for the 
three clouds with adequate spatial resolution to measure 
the diameters. Also included arc the estimated masses of 
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Table 3: CO-H2 Conversion Factors and Cloud Masses 



GMC 


a/a ga i 


M VT 
10 5 Mq 


M ff2 
10 5 Mq 


1+2 
3 

4 
5 


7.2 ± 3.6 
5.9 ± 3.5 
6.7 ± 3.5 


4.8 ± 1.8 
0.87 ± 0.44 
2.1 ± 0.85 


10.9 ± 3.8 

3.5 ± 1.3 

0.64 ± 0.32 

1.5 ± 0.63 



H2. For clouds 1+2 this is obtained using the average con- 
version factor described below, for clouds 3, 4 and 5 this 
is simply the virial mass reduced by a factor of 1.36 for 
the helium contribution. 

The values for the individual clouds are consistent with 
each other, given the rather large error bars. The scatter 
between the values, a, is 0.6, smaller than the individual 
errors, and smaller than the 20% accuracy of the calibra- 
tion of the CO flux. Combining these two contributions, 
we adopt 1.5 as the error for the average conversion fac- 
tor in NGC 1569 instead of using the statistical scatter of 
the three values, thus a/a ga i = 6.6 ± 1.5. We can com- 
pare this with conversion factors calculated with the same 
method by Wilson ( 1995 ) for other dwarf galaxies of nearly 
the same metallicity as NGC 1569. NGC 6822, with 12 + 
log(0/H) = 8.20, has a/a gai < 2.2 ± 0.8, while IC10, 
with 12 + log(0/H) = 8.16 has a/a gal = 2.7 ± 0.5. Thus 
NGC 1569 has a higher conversion rate than dwarf galax- 
ies of similar metallicity by a factor or 2 to 3. 

Israel ( |1997[ ) has used far infrared data from IRAS to 
determine the conversion factor in several magellanic ir- 
regular galaxies, including NGC1569. He obtains a/a ga i 
= 70 + 35 for NGC1569, nearly ten times our value. At 
least two possibilities exist to explain this discrepancy. 
Our value is only valid for those GMCs in which we de- 
tect CO emission, while there may exist regions contain- 
ing molecular material, but no CO. Thus Israel's value 
may be an average over the whole galaxy, including areas 
where there is no CO emission. Alternately, one or more 
of the assumptions used by Israel in his determination of 
the conversion factor may not be valid. For example, his 
method assumes a constant dust-to-gas ratio everywhere 
in the galaxy, which may not be the case. Too little is cur- 
rently known about the distribution of cool dust in dwarf 
galaxies. 

Using our derived conversion factor, the total mass of 
H 2 in the five detected GMCs is (16.5 ± 4.1) x 10 5 M , 
compared to 1.1 x 10 8 M in HI (Reakes [l980| ) . The mass 
of a typical single HI clump that can be identified in the 
map of Israel & van Driel (1990) is of order 10 6 M Q , i.e. 
similar to the mass of the GMCs we find. For comparison, 
we can use the absolute magnitude (Mb = -16.9; Tully 
1988) and an assumed stellar mass-to-light ratio of 1 to 
estimate the total stellar mass to be ~ 10 9 Mq. This mass- 
to-light ratio may be an underestimate, because NGC 1569 
has had a recent burst of star formation, but we see that 



the stellar mass is much larger than the contribution made 
by the molecular or atomic ISM. 

Greve et al. (1996) have not found any CO emission 
in other regions in NGC 1569, but even if there was some 
they missed, it is unlikely to account for much more than 
what is already seen. We conclude that the molecular 
gas contributes a small fraction by mass to the ISM in 
NGC 1569. However, it should be noted that the conver- 
sion factor we derive should only be applied in a strict 
sense to the GMCs. The diffuse emission not detected by 
our interferometric observations may have a different con- 
version factor, since it will have a different (lower) den- 
sity, and possibly a different temperature. 

We can obtain an order of magnitude estimate of the 
density of this diffuse gas. Greve et al. argued that the 
total molecular mass they detected is of the order 2 x 10 6 
Mq for an assumed conversion ratio of ~ 4 times the 
Galactic value. This should give a lower limit to the den- 
sity, because for diffuse gas self shielding of CO is less 
efficient and CO will be more easily dissociated. Thus a 
given mass of H2 will have a higher conversion factor. The 
interferometer detects 22% of the 1— »-0 found by Greve et 
al., leaving a mass of 1.6 x 10 6 Mq undetected. For a 
spherical volume of 120 pc radius (Greve et al.), this cor- 
responds to an average density of ~ 10 molecules cm~ 3 . 

4. Discussion 

4-1. The Physical Conditions of the Molecular Gas 

Line ratios are often used as indicators of the physical 
conditions in the molecular ISM. We will compare the 
12 CO 2— >1 / 1— >0 line ratios we present in Table 2 with 
those found for molec ular gas in various environments. 
Sakamoto et al. ( 1994 ) have derived this line ratio for the 
GMCs Orion A and B, obtaining 0.77 and 0.66 respec- 
tively, similar to our values for GMCs 3 and 4. Thornlcy 



& Wilson ( |1994D have obtained a line ratio of 0.67 ± 0.19, 
averaged over several GMCs in M33. Line ratios signifi- 
cantly higher than unity are often considered indicative 
of warm, optically thin gas and have been observed in in- 
teracting galaxies, or in galaxies with nuclear starbursts 
(Braine & Combes |1992|; Aalto et al. |1995|). For the star- 



burst system M82, Wild et al. (|1992|) find a line ratio of 
1.3 ± 0.3. In NGC1569 the GMCs with the highest line 
ratios (GMCs 1+2) are not significantly nearer to either 
the SSCs or the closest HII regions in projected separation 
than the other clouds, so it is unlikely that warming due 
to star formation is responsible. More likely the line ratios 
reflect a contrast in the densities of the GMCs. 

Because we have data only in the 12 CO 1— >0 and 2— *1 
transitions, we cannot provide tight constraints on the 
physical conditions in the CO emitting gas. Generally a 
line known to be optically thin is necessary for this, such 
as 13 CO 1— >0. Given the low emission in the lines of the 
much more common 12 CO, obtaining 13 CO detections is 
difficult in dwarf galaxies. 
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Fig. 7. Results of an LVG model with 12 CO/H 2 
2 x 10~ 5 and a velocity gradient of 1 km s _1 pc -1 . 



However, we can run a series of large velocity gradient 
(LVG) models and see what ranges of parameter space 
are consistent with our observed 2— >0/1^0 line ratios. 
As input parameters into the models we have a 12 CO/H2 
abundance ratio of 2 x 10~ 5 and velocity gradients of 
0.2, 0.4, 1.0 and 2.0 km s -1 pc -1 . The abundance ratio 
is simply the standard assumption scaled by 0.2, as the 
metallicity of NGC1569 is approximately 20% solar. 

Because we do not have a third line, such as 13 CO 
1— >0, the kinetic temperature Tki n and the density, nn 2 
are degenerate and a large number of solutions may give 
the same line ratios by trading off temperature versus den- 
sity. Figure 7 shows the result of one LVG model, where 
the velocity gradient is 1.0 km s _1 pc _1 . For this set of 
parameters, we see that the gas with the high line ratio 
will have a narrow set of acceptable values for kinetic tem- 
perature (~ 150 K), and densities (log hh 2 ~ 3.5). For the 
low line ratio gas, a wide range of acceptable values exists. 

The physical conditions in the molecular gas deter- 
mine the value of the CO-H2 conversion factor. Because 
the 1— >0 line is optically thick under the conditions found 
in GMCs, the empirical relationship that permits it to be 
used as a tracer of H2 mass depends on a very clumpy 
molecular medium with a low filling factor, where the 
clumps do not shadow one another. If this is not the case 
in NGC 1569 and there is substantial shadowing, then this 
would reduce the derived M mo i, which would, in turn, in- 
crease the conversion factor. 

Alternately, if the GMCs are in a stronger UV radi- 
ation field than those in either IC10 or NGC 6822, the 
CO emission could be reduced due to increased photo- 
dissociation of the CO molecule relative to those two 
galaxies. This would also increase the conversion factor. 
This explanation is consistent with the different natures 



of these three galaxies. NGC 1569 is a BCD which has 
had a major burst of star formation in the recent past, 
which resulted in the two SSCs. The presence of HII re- 
gions also indicates that some star formation is currently 
proceeding. Both IC10 and NGC 6822 are far more qui- 
escent than NGC 1569, although a high concentration of 
Wolf-Rayet stars in IC10 does indicate a recent star forma- 



tion episode in that system (Massey & Armandroff 1995) 



If the increased photoionization is the correct explanation, 
then we would expect to see enhanced emission in the far 
infrared and submillimeter lines of CI and [Gil] compared 
to NGC 6822 and IC10. We will discuss what constraints 
upon the physical conditions in the GMCs are imposed by 
the observed line ratios in the next section. 



4-2. The Relationship Between the Super Star Clusters 
and the GMCs 

For comparison of the CO emission to the optical com- 
ponent of NGC1569, we obtained an HST WFPC2 image 
from the data archive of the Space Telescope European 
Coordinating Facility. Figure 8 shows the CO 1— >0 con- 
tours superposed on this image. The F555W filter used 
corresponds approximately to the V band. Two circles 
show the positions of the SSCs, and crosses show the posi- 



tions of HII regions identified by Waller fll99ip . SSC A is 
near the galaxy center, and about 115 pc east of the near- 
est GMC, number 3. SSC B lies about 53 pc southeast of 
SSC A. No molecular gas is directly associated with these 
clusters. This is to be expected, as these clusters represent 
very strong star formation episodes approximately 15 Myr 
ago. Energy input into the ISM from the by-products of 
star formation (stellar winds and supernovae) may have 
disrupted the natal clouds responsible for the formation 
of the SSCs. Indeed, the HI hole, the extended Ha emis- 
sion, and the hot X-ray gas all attest to the influence the 
SSCs have had upon the ISM. 

An HII region does fall partially within the contours 
of GMC 3, although it is impossible to tell from our data 
if the two are physically associated. In more massive disk 
galaxies, the scale heights of the cold ISM and the young 
stellar populations are small enough that a spatial over- 
lap such as is seen here would be sufficient to assume an 
association between the molecular gas and the HII region. 
However, due to their shallower gravitational potentials, 
dwarf galaxies often have thicker disks than do spirals (e.g. 
Holmberg II, Puche et al. |l~992;). Still, even if NGC 1569 



has an HI scale height of ~ 600 pc like Ho II, the height of 
the molecular gas must necessarily be smaller. Unless they 
are in a non-equilibrium state kinematically, the ensemble 
of GMCs must lie in the plane of the galaxy. The veloc- 
ity dispersion between the clouds (along our line of sight, 
of course) is ~ 7 kms -1 , compared to a global Vf w h m of 



72 kms" 1 for the HI (Reakes |1980p . The scale height will 
be appro ximat ely proportional to the velocity dispersion 
(Kellman 1972 ), which implies a scale height of ~ 60 pc 
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Fig. 8. The contours of 12 CO 1—^0 from Figure 3 super- 
posed over an HST WFPC2 image of NGC 1569 taken 
through the F555W filter. The two circles show the posi- 
tions of SSCs A and B, and the c rosses indicate the posi- 
tions of HII regions from Waller ( 1991 ). 



for the molecular material. With such a low scale height, 
it is likely that GMC 3 and the HII region are physically 
close to each other. The spatial relationship between this 
HII region and the molecular clouds is similar to what is 
seen in the molec ular g as south of 30 Doradus in the LMC 
(Johansson et al. 1998 ), where an HII region lies along the 
edge of the molecular gas, partially overlapping. 

Prada et al. 1994 have identified a star cluster in this 



region of Ha emission and suggested that it might be a 
SSC in the process of forming. Based upon our observa- 
tions, we now consider this to be unlikely. The average 
gas density in the region near this cluster is ~ 200 cm~ 3 , 
assuming a spherical geometry. Sternberg ( 1998] ) has ar- 
gued that the density of molecular gas in the cloud from 
which the SSCs formed was ~ 10 5 cm~ 3 , three orders 
of magnitude higher than what we infer. To recover such 
a high density through clumping of the molecular ISM 
would require an unreasonably large volume filling factor 
of - 0.001. We note that 200 cm" 3 is only a factor of 30 
larger than our lower limit for the density of the diffuse 
molecular gas. 

4.3. The Hot and Cold Phases of the ISM 

NGC 1569 is often cited as a case of a dwarf galaxy ex- 
periencing a blowout of the ISM due to the effects of a 
star formation burst. Ha emission has been found to form 



Fig. 9. The contours of 12 CO 1— s-0 from Figure 3 super- 
posed over an Ha image from Devost et al. 1997 

a halo of emission around the galaxy, with shell struc- 
tures discernable (Hunter et al. 1993 ; Devost et al. 1997 ). 
The radial velocities of as much as ± 200 kms" 1 relative 
to the systemic velocity suggest expanding superbubbles, 
and X-ray data find hot (10 7 K) gas in the interior of these 
bubbles (Heckman et al. |1995| , Delia Ceca et al. |l996| ). The 
mass of the hot X-ray gas is 1.2 x 10 6 f 5 Mg, where f is 
the filling factor of the gas. It is reasonable to expect that 
the violent process of heating the gas and driving it in 
an outflow would leave some kind of observable signature 
upon the remaining cold ISM. 

The HI hole centered on SSC A discovered by Israel 
& van Driel ( 1990D could be an example of this. They ar- 
gue that the data are consistent with a picture in which 
the formation of the hole began about 10 7 yr ago, driven 

by the expansion of supernova remnants. The angular size 

// 
they measure for the hole is ~ 10. , with the result that 

the GMCs we have imaged lie just outside the hole. The 

expanding Ha bubbles have a dynamical age of ~ 10 yr 



(Heckman et al. 1995), similar to the age of the HI hole. 



This provides a limit to the duration of the starburst that 
created SSC A. Based upon the shape of the n on-th ermal 
radio continuum spectrum, Israel & de Bruyn ( 1988 ) have 
argued that the star formation burst in NGC 1569 ended 
about 5 x 10 6 yr ago, which would then suggest a burst 
duration of about 5 x 10 6 . There are three possible sce- 
narios regarding a connection between the hot and cold 
phases of the ISM. 

1 . The GMCs are a direct result of the expansion of the 
hole in the HI driven by the hot gas. A shell of ac- 
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cumulated material along the edge of the hole might 
form molecular clouds. Indeed, an expanding shell of 
CO emission has been detected around the giant HII 
region 30 Doradus by Cohen et al. ( |1988| ). But in the 
case of NGC 1569, the GMCs do not form a shell-like 
distribution around the HI hole, nor do the kinematics 
of the CO emission indicate an expanding shell, so this 
scenario is not very likely. 

Instead of originating in a swept up shell, the GMCs 
may have collapsed from pre-existing high density ma- 
terial due to the shock of the outflow. The HI maps of 
Israel & van Driel do show that the peak of the HI col- 
umn density is west of the SSCs. The lack of molecular 
gas elsewhere in the vicinity of the SSCs (Greve et al. 



1996J) would be explained if the HI elsewhere did not 



have a high enough density for GMCs to collapse even 
with the catalyst of a passing shock. In this case, the 
CMC formation and subsequent star formation would 
be a prime example of star formation being triggered 
by an earlier, nearby star formation event. 
3. There is no relationship. The GMCs could predate the 
star formation burst of about 15 Myr ago that created 
the SSCs. Molecular clouds will tend to form where 
the gas density is high. A number of dwarf irregular 
galaxies without large star formation bursts are known 
to have irregular, clumpy HI distributions (e.g. Sag 
DIG, Young & Lo [l997|; GR8, Carignan et al pJ90| ), 
and near such clumps is the natural place to expect 
molecular gas. 

Option 1 is unlikely, for the reasons stated above. Op- 
tion 2 is more likely, given the relative geometry of the 
SSCs and the GMCs. An interesting morphological note 
is the Ha arm identified by Waller ( J1991 ). This feature 
extends about 640 pc from the main part of NGC 1569 
and it connects to the galaxy very close to the position 
of the GMCs and cluster C. A similarly shaped feature 
occurs in HI, sitting just to the exterior of the Ha feature. 
Waller interprets these as the interface between outflowing 
hot gas and the cool neutral material. The GMCs would 
then be positioned at the part of this interface region with 
the highest gas density. This coincidence makes us favor 
option 2, although the evidence is certainly not conclu- 
sive. Further CO observations of otherwise similar dwarfs 
which lack SSCs will help us to understand this issue. We 
discuss high resolution CO observations of other BCDs in 
section 4.4. 

If we assume a star formation efficiency (SFE), de- 
fined as the fraction of gas mass converted into stars, we 
can estimate the amount of gas required in the burst that 
created the SSCs. Sage et al. ( |1992| ) have calculated global 
SFEs for a number of BCDs using Ha, CO and HI obser- 
vations. These SFEs are calculated as the star formation 
rate (from Ha observations) per gas mass - i.e. the in- 
verse of the gas consumption timescale. To arrive at our 
definition of SFE, we must multiply by a burst duration, 



which we will take to be 5 x 10 6 yr. Note that implicit 
in this is the further assumption that the current star for- 
mation rate is the same as the average star formation rate 
during the burst. We will take an average over several 
BCDs, with the assumption that since they will be at dif- 
ferent stages in the development of their star formation 
episodes, the result will be approximately an average SFE 
for the duration of a typical burst. Averaging SFEs for the 
galaxies which most resemble NGC 1569 (which we define 
as having an HI mass within a factor of 2 of NGC 1569), 
we obtain 1.9%. If the burst duration is longer than we 
have assumed, then the SFE will increase, because a larger 
fraction of the gas mass will have been converted to stars 
using the average star formation rate. For a duration of 
~ 10 8 yr, the SFE will approach 100%, much higher than 
is expected. However, burst durations much longer than 
5 x 10 6 yr confl ict wi th the dynamical age of the HI hole 
(Heckman et al. 1995|) and the time for th e end of the star 
formation burst (Israel & de Bruyn 



The SFE derived above compares well with the values 
of 2.4% and 1.9% determined by counting individual stars 
in the giant HII regions NGC 595 and NGC 604 in M33 
by Wilson & Matthews ( 1995[ ). Of course these two are 
individual giant HII regions, and thus their SFEs are not 
global values, as come from Sage et al. But the similarity in 
the SFEs between the two different environment obtained 
using two different methods suggests that the values are 
reasonable. 



De Marchi et al. (1997) have estimated the mass of 
SSC A at 2.8 x 10 5 M , so a SFE of 1.9% yields an 
original gas mass (H 2 + HI) of 1.5 x 10 7 Mq. Wilson 
& Matthews (1995) find the ratio of molecular to atomic 
hydrogen in NGC 595 and NGC 604 to be approximately 
1:1, so this would imply an original Mh 2 of 7.5 x 10 6 
Mq for the gas that formed SSC A. This is larger than 
the 1.7 x 10 6 M Q found in the current GMCs. Of course 
there is significant diffuse emission that was not detected 
in our interfero meter observations. In the center pointing 
of Greve et al. ([1996 ) this amounts to nearly a factor of 
4.5, so we can estimate a lower limit on the total H 2 mass 
in our field to be 7.7 x 10 6 . This is similar to the estimate 
for the gas that created SSC A. Thus there is sufficient 
gas present now to explain the SSCs, but distributed over 
an area ~ 200 pc in diameter. This extended distribution 
of the gas may explain why we see current star formation 
as typical HII regions, but not as newly born SSCs. 

4-4- Comparision with Other Galaxies 

AAA. dlrrs 

Several of the irregular galaxies in the Local Group have 
been observed with high spatial resolution in CO, includ- 
ing the LMC and SMC (each observed with SEST), and 
IC10 and NGC 6822 (both observed at OVRO). When in- 
dividual molecular clouds are resolved, they follow a size- 



10 



C.L. Taylor et al.: Giant Molecular Clouds in the Dwarf Galaxy NGC1569 



linewidth relationship very similar to that of Milky Way 
GMCs. The clouds seen in these more nearby galaxies tend 
to be somewhat smaller that those in NGC 1569, e.g. ~ 30 



pc diameter in the LMC and SMC (Johansson et al. 1998; 
Rubio et al. 1993). It is likely that higher resolution obser- 
vations would separate the NGC 1569 clouds into smaller 
units, as it is known that molecular clouds are clumpy and 
have a low volume filling factor. The number of clumps 
contained in any arbitrary structure is not as important 
as whether or not that structure is gravitationally bound, 
and in virial equilibrium. 

4.4.2. BCDs 

Except for a few cases, the history of observing CO in 
BCDs is largely one of non-detections (e.g. Young et al. 



et al. 1998 and Gondhalekar et al. 1998). Because BCDs 



1986| , Israe l fc B urton [1986, Tac coni & You ng |1987| , Ar- 
nault et al |1988| , Sage et al. 1992| , Is rael e t al. |1995| , Taylor 



have comparatively high star formation rates, the lack of 
detections is not likely caused by a lack of molecular gas, 
which has been a frequent, but false, conclusion in the 
past. Instead, it is more likely to be attributed to the gen- 
erally low metallicities of most BCDs (Searle & Sargent 



1972) leading to a high CO-H2 conversion factor. 



Even fewer BCDs have been mapped in transitions 
of molecular gas. These include NGC4214 (Becker et al. 



1995), NGC5253 (Turner et al. 1997), Henize 2-10 (Kob 



ulnicky ct al. [1995[ , Baas et al. |1994[), Mrkl90 (Li et al. 

We shall discuss 



1994) and III Zw 102 (Li et al. 



1993b) 



these in the remainder of this section to put our results 
on NGC 1569 into a larger perspective. 

NGC4214: Becker et al. ( |1995| ) have mapped NGC4214 
simultaneously in the 1— >0 and 2^0 lines of 12 CO using 
the 30-m telescope of IRAM. They detected a large region 
of emission, about 1000 pc x 700 pc in size, near the cen- 
ter of the galaxy. This emission shows structure on scales 
of ~ 500 pc, which was approximately the resolution limit 
of those observations. One feature which is well resolved 
has a virial mass of ~ 10 7 Mq, larger than what we have 
seen in NGC1569. Clouds of this size are rare in surveys of 
Galactic GMCs (e.g. Sodroski |l99l| ; Sanders et al. |l985[ ), 
so perhaps these features in NGC4214 are simply collec- 
tions of unresolved smaller clouds. When observed with 
a large enough beam, the GMCs we see in NGC1569 do 
appear as a si ngle l arge cloud of diameter ~ 150 - 200 pc 
(Greve et al. (|1996|) . 

Becker et al. find a 2^0/1^0 line ratio of 0.4 ± 0.1 
for NGC4214. Their data have relatively low spatial res- 
olution (13" in 12 CO 2->l) and NGC4214 is more than 
twice as far away as NGC1569. Therefore their line ratio 
represents not a value for an individual cloud, but an av- 
erage over multiple clouds belonging to a molecular cloud 
complex. In addition, they used a single-dish telescope, so 
they do not have the problem of missing flux due to a lack 
of short spacings. It would be interesting to obtain inter- 



ferometric CO observations of NGC4214 in order to derive 
line ratios on smaller physical scales than was possible for 
Becker et al. and see if any dense clouds with high line 
ratios are present, such as we find in the case of GMC 
1+2 in NGC1569. 



NGC5253: Turner et al. (|1997| ) have mapped NGC5253 
with OVRO at resolution of 190 x 90 pc (for the distance 
of 4.1 Mpc). Individual GMCs thus are not resolved in this 
galaxy. The CO distribution is only marginally resolved at 
best, and is weakly detected. Turner et al. recover approx- 
imately one half the flux detected by single-dish observa- 
tions (e.g. Taylor et al. 1998). The CO emission is found 
near the optical center of the galaxy, but directly above 
it. It also lies perpendicular to the optical major axis, and 
along a dust lane. Turner et al. suggest that the CO may 
have been ccreted onto NGC5253 from another system. 

Henize 2-10: At a distance of approximately 9 Mpc, He 
2-10 is too distant to resolve even molecular complexes, as 
was done for NGC4214. However, global line ratios have 
been obtained by Baas et al. ( 1994 ), who find a 2— >1/1— >0 
line ratio of 0.97 ± 0.16. Because of the low spatial reso- 
lution of their observations and the large distance of the 
galaxy, this likely represents an average over different re- 
gions of the galaxy, with gas in different physical condi- 
tions. Indeed, Baas et al. explain this line ratio, as well as 
the 3-^2/2— +1 line ratio of 1.34 ± 0.17, as resulting from 
a two temperature model with a component of the CO 
emitting gas at a temperature of < 10 K, and another at 
> 75 K. Because we only have data in two transitions, we 
cannot constrain sophisticated models, but we certainly 
cannot exclude such a two temperature model for the CO 
emission in NGC1569. 

The spatial distribution of 12 CO 1— >0 emission in He 2- 
10 is described by Kobulnicky et al. (1995), who obtained 



interferometer observations with OVRO. The peak of the 
CO emission is located a few arcseconds from the regions 
of current star formation, which is consistent with what 
we find in NGC1569. However, He 2-10 has an unusually 
extended CO distribution, with a spur of CO emitting gas 
extending southeast from the star forming center of the 
galaxy. This feature is also reproduced in the HI data. 
Kobulnicky et al. suggest that He 2-10 is a moderately 
advanced merger between two dwarf galaxies. He 2-10 has 
a hole in the HI distribution near the current star-forming 
regions, but unlike NGC1569, this hole is filled by the bulk 
of the observed molecular gas, and there is no evidence for 
an outflow of hot gas from the starburst region. 

Mrkl90: This galaxy was observed with the OVRO in- 
terferometer in the 12 CO l->0 line by Li et al. (|1994|). Al- 



though the CO emission is only marginally resolved (dis- 
tance = 17.0 Mpc), Li et al. find evidence that the molec- 
ular gas is distributed in a ring centered on the galaxy 
center. They suggest that starbursts in the central region 
of the galaxy have acted to clear gas out from this area, in 
much the same way as is seen in NGC1569. There is, how- 
ever, no direct evidence for this process in Mrkl90, unlike 
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NGC1569. Li et al (1993a) report single-dish observations 
in multiple transitions of CO for this galaxy. They obtain 
a 2— >1/1— >0 line ratio of 0.93 ± 0.25, which is roughly 
consistent with the line ratios for NGC1569. 

Ill Zw 102: Despite its relatively high optical luminos- 
ity (Mb = -19.2), this galaxy is often included in studies 
of dwarf galaxies (e.g. Thuan & Martin |1981| ). Li et al. 
1993b| obtained 12 CO l->0 maps with OVRO, and also 
single dish spectra in the 2— >1 and 1— >0 transitions with 
the IRAM 30-m telescope. The distance to III Zw 102 is 
large (23.5 Mpc), so structures like GMCs could not be 
resolved in the observations. They found a 2— >1/1— >0 line 
ratio of 0.66 ± 0.12, quite similar to what we found for 
GMCs 3 and 4 although, given the large errors on our line 
ratios, the value for GMC 1+2 is also consistent. This line 
ratio for III Zw 102 is measured over the 23" beam of 
the 30-m telescope at 115 GHz. This corresponds to ~ 2.6 
kpc within the galaxy, so the line ratio is clearly an aver- 
age over a large number of GMCs. The distribution of the 
CO emitting gas is roughly correlated with the star for- 
mation as traced by Ha and radio continuum images, but 
the peak of the CO emission is offset relative to the peak 
of the optical emission. This is quite similar to NGC1569, 
in which we find molecular material associated with HII 
regions, but not with the SSCs. 



5. Summary and Conclusions 

We have presented CO 1— >0 and 2— >1 observations of 
the dwarf starburst galaxy NGC1569 obtained with the 
IRAM millimeter interferometer. We confirm the result of 



Greve et al. (1996) that the molecular gas is not associ- 
ated with the super star clusters in NGC1569, but instead 
with nearby HII regions. The major results of our study 
are: 

1. The CO emission is resolved into a number of individ- 
ual giant molecular clouds. These GMCs have sizes and 
linewidths similar to the more massive GMCs in the 
Milky Way Galaxy, and to those in other Local Group 
galaxies. The 2— >1/1— >0 line ratios of the GMCs in 
NGC1569 range from 0.64 ± 0.30 to 1.31 ± 0.60. The 
lower values are similar to what is typically seen in the 
Milky Way GMCs, while the higher values are similar 
to what is observed in starburst galaxies. 

2. The CO-H2 conversion factor is found to be 6.6 ± 1.5 
times the Galactic value by applying the virial theo- 
rem to 3 of the GMCs detected in NGC 1569. This 
is approximately three times higher than is found for 
NGC6822 and IC 10, two dwarf irregular galaxies in 
the Local Group with nearly the same metallicity as 
NGC 1569. This difference may due to a stronger UV 
radiation field in NGC1569 compared to the dwarf ir- 
regulars. Sensitive observations of cooling lines like CI 
and CII, expected to be enhanced in the presence of 
strong photo-dissociation, are called for to decide this 



matter. NGC1569 has recently experienced a strong 
burst of star formation which formed the two super 
star clusters, and still has a number of HII regions. 
3. The GMCs are observed to be just outside the edge of 
the HI hole surrounding SSC A. This hole is thought 
to be a region swept clear of cold gas by the observed 
outflow of hot X-ray emitting gas. It is possible that 
shocks from this process may have contributed to the 
formation of the GMCs, although our data do not place 
any constraints on this scenario. 
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